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vi è quella di permettere di stabilire quando una matrice è invertibile, cioè 
quando, data una matrice A , esiste una matrice B tale che AB = id = BA. 

La definizione stessa di determinante e la dimostrazione delle molte, rile- 
vanti proprietà del determinante sono tecniche e abbastanza noiose. Allora 
daremo, come definizione di determinante, un algoritmo per calcolarlo e un 
elenco delle proprietà assolutamente da non ignorare. Raggrupperemo in una 
parte successiva le dimostrazioni. 

Definizione 3.5.1 Sia A una matrice quadrata ( N.B. Il determinante è 
definito solo per matrici quadrate). 

1) Se V ordine di A è 1, cioè A è costituita da un solo elemento A = (a), 
a E M, allora , per definizione, dei A = a. 

2) Se A ha ordine n > 1, allora dei A si calcola secondo la seguente formula 
ricorsiva : 


det A 


ctn det Ah ~ &12 det Ai 2 + • • • + (— det Ai & + • • • + 
+( — l) 1+n &in det Ai n , (3.5) 


dove con Aij si indica la matrice quadrata di ordine n — 1 ottenuta da A can- 
cellando la i-esima riga e la j -esima colonna, detta complemento algebrico 
deir elemento a^. 


La formula utilizzata per questa definizione viene detto sviluppo di Laplace 
del determinante secondo la prima riga. 

a b 
c d 


Sia A di ordine 2, A = 


Applicando la definizione: 


det A = a det Ah — b det A 12 . 


Ma An = (d) e A 12 = (c) e per il punto 1) della definizione si ha: 

det A = ad — bc . 


Esempio 3.5.2 A = ^ ^ det A = 1 — (—2) = 3. 

Ormai risulta chiaro che, se ho una matrice di ordine n, sviluppando secondo 
Laplace, cioè secondo la formula (3.5), ho una somma con n determinanti di 
matrici di ordine n — 1 che, a loro volta, vanno sviluppati allo stesso modo, 
per un totale di n{n — 1) determinanti di matrici di ordine n — 2 e così via 
fino ad arrivare a dover calcolare ni determinanti di matrici di ordine 1 per 
le quali si applica la prima parte della definizione del determinante. 
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(3 1 2\ 

Esempio 3.5.3 A = 0 —1 2 I . 

Vi i i j 

ì 1 f) + (~l) 1+2 det (° ^ +(-l)‘+ 3 2det 
= 3(-3) - (-2) + 2 • 1 = -5. 

Principali proprietà del determinante 

1) Lo sviluppo di Laplace del determinante può essere fatto secondo una 
qualsiasi riga o colonna, non solo la prima. Se scelgo la i-esima riga, la 
formula diventa: 

det Li = (— l)* +1 aji det A a + (— l) ,+2 aj 2 det A i2 + • • • + 

+ ( — 1 Y +k aik det Aik + ■ • • + (— 1 ) t+n cii n det Ai n 

n 

= y^(-l y +k aik det A ik ■ 

k = 1 

Se scegliessi la j-esima colonna, farei variare da 1 a n l’indice di riga e terrei 
fisso j come indice di colonna: 

n 

det A = det A k j . 

k = 1 

Il risultato non dipende dalla scelta della riga o della colonna rispetto a cui 
faccio lo sviluppo. 

2) Se una matrice ha una riga (o una colonna) fatta tutta di zeri, il suo 
determinante è zero. 

3) Se una matrice A' è ottenuta scambiando tra loro di posto due righe (o 
due colonne) della matrice A , allora det A' = — det A. 




4) Se una matrice ha due righe (o due colonne) uguali, non necessariamente 
adiacenti, allora il suo determinante è zero. 
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5) Se A! è una matrice che si ottiene moltiplicando per uno scalare À tutti 
gli elementi di una riga (o di una colonna) della matrice A, allora: 

det A' = X det A . 

6) Se À è uno scalare, det (ÀA) = À n det A, dove n è l’ordine di A. 

7) Siano Ri e i? 2 due righe (risp. due colonne), cioè due matrici di tipo 
lxn (risp. nxl). Sia A una matrice di ordine n che abbia come i-esima riga 
la riga Ri + R 2 e siano Ai e A 2 due matrici di ordine n aventi come i-esima 
riga rispettivamente Ri e i? 2 e tutte le altre righe uguali alle righe di A e 
nella medesima posizione. Allora: 

det A = det A 1 + det A 2 . 

8) Se si somma ad una riga della matrice A (risp. ad una colonna) una 
combinazione lineare delle altre righe (risp. colonne) di A, il determinante 
della matrice così ottenuta è uguale a det A. 

Ricordiamo che ad esempio le colonne sono matrici di tipo nxle quindi 
di esse si possono considerare le combinazioni lineari, secondo la Definizio- 
ne 3.2.6. 

9) Se I n è la matrice identità di ordine n, det I n = 1, indipendentemente 
da n. 

10) Il determinate di una matrice è uguale a quello della sua trasposta: 

det ( A) = det A . 

Queste proprietà vanno tenute sempre presenti, perché possono facilitare 
enormemente il calcolo dei determinanti. Valga come esempio, il suggeri- 
mento di sviluppare il determinante secondo la riga o la colonna che contiene 
il maggior numero di zeri, poiché il risultato, in base alla proprietà 1) è lo 
stesso, ma ogni zero risparmia il calcolo di un determinante di ordine n— 1. 

Osserviamo che nessuna indicazione viene data per la somma di matrici; 
tra det (A + £>), det A e det B non vi è nessuna relazione significativamente 
utilizzabile. Vi è invece una importante relazione per il determinante del 
prodotto tra matrici. 
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Teorema 3.5.4 (Binet) Il determinante del prodotto di due matrici è ugua- 
le al prodotto dei determinanti dei fattori; in breve: 

det (A • B) = det A det B . (3.6) 

Esistono molte dimostrazioni di questo teorema, spesso tecniche e poco istrut- 
tive. Ne daremo daremo una un po’ articolata ma che ci sembra utile. 
Premettiamo la seguente: 

Osservazione 3.5.5 II determinante di una matrice triangolare si calcola 
facendo il prodotto degli elementi sulla diagonale principale. Basta infatti 
sviluppare il determinante secondo la prima colonna e ripetere il procedimento 
per il complemento algebrico che è ancora triangolare. 

S e A e B sono triangolari, per TOsservazione 3.3.8 sul prodotto di matrici 
triangolari, Ti-esimo elemento sulla diagonale principale del prodotto C è 
uguale al prodotto degli i-esimi elementi delle diagonali principali di A e B 
rispettivamente: Cu = a^ò^. Quindi, usando TOsservazione 3.5.5, abbiamo 
che: 

det (A - B) = a n b u * • • a nn b nn = a n • • • a nn b u • • • b nn = det A det B . 

Nella sezione 3.4 (Proposizione 3.4.8) abbiamo visto che è possibile trasfor- 
mare una matrice in una triangolare tramite trasformazioni elementari di riga 
del secondo tipo. Si noti che, per la proprietà 8) del determinante, queste 
trasformazioni non alterano il valore del determinante (a differenza di quelle 
del primo tipo, che lo cambiano di segno, e di quelle del terzo tipo, che hanno 
il risultato di moltiplicarlo per uno scalare). 

Adepereremo questo risultato per dimostrare il teorema di Binet nel caso 
di matrice generiche (non necessariamente triangolari). 

Dim. (del Teorema di Binet). Date A e £>, siano T R (A) e Tc(B ) le matrici 
equivalenti per righe e risp. per colonne ad A e B che siano triangolari e ot- 
tenute con trasformazioni del secondo tipo, quindi con lo stesso determinante 
di A e £>, rispettivamente. Allora: 

det (AB) = det (T C (T R (A • B))) = det(T c (T R (A) •£)) = det(T R (A) -T C (B)) . 

Questo applicando per la prima uguaglianza le proprietà del determinante e 
delle trasformazioni elementari, e per le altre la Proposizione 3.4.7. Poiché 
T r (A) e Tc(B ) sono triangolari, per TOsservazione 3.3.8, 

det (T r (A) ■ T C (B)) = det (T R (A)) det (T C (B)) . 
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Infine, poiché det (T R (A)) = det A e det (Tc(B)) = det £>, si ha la tesi: 
det (AB) = det A det B . C 

Dim. (delle 10 proprietà del determinante). 

1) Dimostriamolo per induzione sull’ordine della matrice A. Se A ha ordine 
1, non vi è scelta, la proprietà vale. Supponiamo che valga per ogni intero 
minore di n e mostriamo che vale per n. Lo sviluppo del determinante secondo 
la i-esima riga è dato da: 

(— l) 2+1 a^i det Ah + • • • + (— l) z+fc a^ det A^ + • • • + (— 1 ) t+n a m det A in . 

A questo punto sviluppo il determinante di ciascun A^ secondo la prima riga 
e ottengo una sommatoria su h e k di termini del tipo: 

(-1 T +k+1+h digiti det (Afe) ih 

con i^leh^ke dove (Afa) ih è la matrice che si ottiene da A cancellando 
la prima e la z-esima riga e la h-e sima e la /c-esima colonna. 

Se sviluppiamo il det A secondo la definizione, cioè la prima riga, si ha: 

dii det An + * * * + (—l) 1+h di h det A\h + • • • + (— l) 1+n ai n det A\ n . 

Se sviluppiamo il determinante di ciascun Aih secondo la sua (z — l)-esima 
riga che corrisponde alla i-esima riga di A , per ipotesi induttiva il risultato 
non cambia, ma formalmente otteniamo una somma di termini del tipo: 

(-1 T +k+1+h dikdih. det (AAfe 

con i^leh^ke dove (Aih)ik è la matrice che si ottiene da A cancellando 
la prima e la z-esima riga e la à-esima e la /c-esima colonna e dove i termini 
ciik sono riferiti ad A , nel senso che il termine che si trova al posto (i — 1, k) 
di Aih è proprio il termine che si trova al posto (z, k) di A, se k < /i, e al 
posto (i, k + 1), se k > h. 

Con analogo ragionamento, cioè prima sviluppando per una colonna e poi 
per una riga, quindi sviluppando per una riga e poi per una colonna, sempre 
per induzione, si dimostra che il calcolo del determinante sviluppato secondo 
una colonna è lo stesso che sviluppato secondo una riga. 

2) Ovvio, a questo punto, perché sviluppando secondo la riga, o colonna, 
tutta di zeri, l’espressione che si ottiene è zero. 

3) Supponiamo che A' sia la matrice tale che la z-esima riga di A' è uguale 
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alla (i + l)-esima riga di A e la (i + l)-esima riga di A! è uguale alla i-esima 
riga di A , cioè A' è la matrice che si ottiene da A scambiando due righe 
adiacenti. 

Sviluppando det A secondo la i-esima riga si ha: 

det A = ( — l) iJrl an det Ah S - * • • ~h ( — i y jJrn d^ n det A^ n . 
Sviluppando det A! secondo la [i + l)-esima riga si ha: 
det A' = (— l) (z+1)+1 a'( i+1)1 det H h (-l) (z+1)+n a( i+1)n det . 

Ma, poiché a' (i+1)fc = a ik e A' (i+1)fc = A ik per ogni k = 1, . . . , n, allora 

det A' = (—1) ((— l) z+1 aa det An + • • • + (— l) lJrn a in det A in ) = — det A . 

Se le due righe non sono adiacenti, la matrice che le ha scambiate si ottiene da 
A mediante una successione di matrici Ai ognuna delle quali si ottiene dalla 
precedente scambiando due righe adiacenti, nel modo seguente. Supponiamo 
che voglia scambiare la riga h con la riga h + k, cioè tra le due righe da 
scambiare ve ne sono k — 1. Allora Ai è la matrice che si ottiene da A 
scambiando la riga h con la riga h + 1, poi scambio la riga h + 1 con la h + 2 
e dopo k — 1 scambi tra righe adiacenti ho che la vecchia riga h è adiacente 
alla riga h + k. Le scambio, cosicché adesso la vecchia riga h + k precede la 
riga /i, e dopo altri k — 1 scambi porto la riga h + k ad occupare la h-e sima 
posizione. In totale la matrice A' è stata ottenuta da A dopo 2(k — 1) + 1 
scambi. Siccome ad ogni scambio tra righe adiacenti il determinante cambia 
segno, si ha che: 

det A' = (— l) 2 ^- 1 det A = — det A . 

Stesso discorso per lo scambio di due colonne. 

4) Se A ha due colonne uguali e se A' è la matrice che ha le due colonne 
scambiate, A' = A. Ma anche det A' = — det A. Quindi det A = det A' = 
— det A, che implica det A = 0. 

5) La dimostrazione è immediata. Basta sviluppare det A' secondo la riga o 
la colonna moltiplicata per À, raccogliere À ed accorgersi che ciò che rimane 
è det A. 

6) Applicando n volte il risultato precedente, si ottiene la relazione cercata. 

7) Sviluppando det A secondo la i-esima riga si ha: 

(— l) z+1 a^i det Ah + • • • + (— l) z+n a^ n det A in . 
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Ma a ik = (di)ik + (a 2 )iÀ:j dove (cLj) ik è il corrispondente elemento di Aj, con 
j = 1, 2. Inoltre è facile convincersi del fatto che A ik = (Aj) ik , dove (Aj) ik è 
il corrispondente complemento algebrico di Aj. Abbiamo dunque: 

det A = 

= (— l) z+1 a^i det Ah + • • • + (— l) 2+n a^ n det A^ n 

— ( — l) z+1 (( a i)u + (^ 2 ) 21 ) det Ah + • • • + (— l) z+n ((&i)m + (& 2 )m) det A{ n 

— (— 1) z+1 (^i)ìi det Ah + • • • + (— l) i+n (^i)m det A^ n + 

+( — 1 )* +I (u 2 )u det A^i + • • • + (— l)^ +n (u 2 )m det A^ n = det Ai + det A 2 . 


8) Se alla i-esima riga di A sommiamo una combinazione lineare delle altre 
righe con coefficienti A*, applicando la 5) e la 7) abbiamo che il determinante 
della matrice A' così ottenuta è uguale al determinante di A più una combi- 
nazione lineare con gli stessi coefficienti A* di determinanti di matrici aventi 
due righe uguali, e quindi tutti nulli. 

Procedimento analogo per le combinazioni lineari di colonne. 

9) Per induzione. 

li = (1), det (1) = 1 . 

Supponiamo che det/ n _i = 1. Sviluppando det I n secondo la prima riga si 
ha che: 

det I n — 1 • det 7 n _i = 1 . 


10) Segue subito dalla 1). 


3.6 Invertibilità di una matrice e metodo di 
Gauss 

A questo punto vogliamo affrontare il seguente problema: 

Sappiamo che rinsieme M n (R) delle matrici quadrate di ordine n a coef- 
ficienti reali è un anello (non commutativo) con identità, con le operazioni di 
somma e di prodotto riga per colonna. E esso un campo? Ossia, data una 
matrice non nulla A, esiste una matrice B tale che sia A • B — B • A — I n l 
Il teorema di Binet ci dà subito la risposta. No, se n > 1. Infatti necessa- 
riamente, se A • B = 7 n , det A det B = det (AB) = det7 n = 1. Quindi, se 
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det A = 0 ciò è impossibile. Ed esistono molte matrici A 7 ^ 0 tali che, però, 
det A = 0. Esempio: A = 

Abbiamo allora che, condizione necessaria affinché una matrice quadrata 
A ammetta la matrice inversa, è che det A 7 ^ 0. In realtà tale condizione è 
anche sufficiente. 



Proposizione 3.6.1 Una matrice A G M n (R) è invertibile, cioè esiste B 
(dello stesso ordine) tale che A - B = B • A = I n , se e solo se det A 0. 


Dim. Il ”solo se” lo abbiamo già visto poco sopra. Dimostreremo che se 
det A yé 0, allora esiste la matrice inversa, indicata d’ora in poi con A -1 , 
costruendola esplicitamente. 

Sia Aij il complemento algebrico dell’elemento di posto i,j di A. Sia B 
la matrice tale che il suo elemento bij sia così definito: 




det Aij 
det A 


Osserviamo che qui interviene l’ipotesi che det A 7 ^ 0. 

Verifichiamo che l B — A -1 , cioè f B • A = A • f B = I n . 

Moltiplichiamo la z-esima riga di f B , quindi la z-esima colonna di B , per 
la j -esima colonna di A e vediamo che tale prodotto è zero se i 7 ^ j e vale 1 
se i — j. 



Q'ijbu + * * * + &kjbkì + * * * + a n jb 


ni 


_ ( — det Au + * * * + (— 1 ) k+l akj det A^ + • • • + (— 1 ) nJrl a n j det A n i 

det A 

Se i = j, il numeratore è lo sviluppo di Laplace di det A secondo la z-esima 
colonna. Pertanto si ha che: 



1 


det A 


det A 



Se z yé ji, l’espressione tra parentesi quadra è lo sviluppo di Laplace del 
determinante di una matrice A', che ha tutte le colonne uguali a quelle di 
A tranne la sua z-esima che coincide non con la z-esima di A bensì con la 
j-esima colonna di A. Pertanto A' è una matrice che ha due colonne uguali, 
la z-esima e la j-esima, entrambe uguali alla z-esima colonna di A, e quindi 
det A = 0. Perciò, se z 7 ^ j, Qj = ^^-detA' = 0. Analogo ragionamento 
posso seguire se moltiplico la z-esima riga di A per la j-esima colonna di f B 
che sarebbe la j-esima riga di B. □ 
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Esempio 3.6.2 Sia A 

bile. 


Z 1 2 A 

0 1 1 , det A = | 0 e quindi A è inverti- 

VI - 1 0/ 


&n 
&12 
b 1 3 

Ò 2 1 
&22 
&23 
b-.i i 
&32 
^33 


(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 

(- 1 ) 


2 det Ah 

det A 

3 det Ai2 

det A 

4 det A13 

det A 

3 det -A 2 i 

det A 

4 det A22 

det A 

5 det A23 

det A 

4 det A31 

det A 

5 det A 32 

det A 

6 det A33 

det A 



1 

3 ’ 


B = ( 6 «) 







5z verifichi che A - A 1 = A 1 • A = I n . 


Osservazione 3.6.3 1) La matrice inversa di A se esiste, è unica. Infatti, 
siano B e C tali che A • B = B • A = I n e A • C = C • A = I n . Allora: 

B = B-I n = B- A- C = I n -C = C. 

2) Se A e B sono invertibili, anche A • B lo è e inoltre: 

(A - B)- 1 = B- 1 - A' 1 . 
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Cioè V inversa del prodotto è uguale al prodotto delle inverse ma in ordine 
inverso. Infatti A • B • B~ l • A -1 = A • I n • A -1 = I n . 

Osserviamo che abbiamo sempre utilizzato la proprietà associativa del pro- 
dotto riga per colonna. 

3) ( t A )~ 1 = t (A~ 1 ). Cioè l'inversa della trasposta è la trasposta dell'inversa. 
Infatti: l A • \A~ l ) = t A~ 1 • A = t I n = I n - 

Si può quindi scrivere senza ambiguità t A 1 . 

Metodo di Gauss per il calcolo della matrice inversa. Il seguente 
metodo fornisce un criterio per stabilire se una matrice quadrata è invertibile 
e, in caso affermativo, permette di calcolare la matrice inversa. 

Data una matrice quadrata A, consideriamo la matrice M di tipo n x 2n 
le cui colonne sono ordinatamente quelle di A e quelle della matrice identità 
di ordine n, cioè M è la matrice ottenuta da A affiancandole la matrice I n . 
Scriviamo: 

M = (A,/ n ). 

Per la Proposizione 3.4.8, possiamo ridurre la matrice A ad una matrice 
triangolare A' con trasformazioni elementari di riga. Si noti che A è invertibile 
se e solo se lo è A'. D’altra parte A' è invertibile, cioè det A' ^ 0, se e solo 
se gli elementi sulla sua diagonale sono tutti non nulli. 

Se questo non avviene abbiamo concluso che A non è invertibile. In caso 
contrario possiamo procedere al calcolo dell’inversa di A. Osserviamo che, 
moltiplicando ogni riga per un multiplo opportuno, possiamo allora supporre 
che tutti gli elementi siano uguali a uno. Possiamo poi, sempre con trasfor- 
mazioni elementari di riga, fare in modo che da avere degli zeri sopra ciascuno 
di essi, cioè ridurci ad avere la matrice identità. Il procedimento è simile a 
quello usato nella dimostrazione della Proposizione 3.4.8 (iniziando dall’ul- 
tima colonna anziché dalla prima). Facendo tutte queste trasformazioni di 
riga alla matrice M (non solo su Al) si ottiene allora una matrice 

M' = (J n , B) . (3.7) 

Diciamo che B è proprio l’inversa di A. Infatti le trasformazioni di riga su 
M corrispondono esattamente a moltiplicarla a sinistra per un’opportuna 
matrice T. Abbiamo cioè: 


M' = T • M = T • (A, I n ) = (T • A, T ) , 
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da cui segue per confronto con la (3.7): 


T ■ A — I 


n > 


T = B 


quindi B • A = I n . 


3.7 Caratteristica di una matrice 


Per matrici di tipo qualunque, anche rettangolari, si può introdurre un altro 
indicatore che risulterà di fondamentale importanza. 

Definizione 3.7.1 Data una matrice A , si dice minore di A (di ordine k) 
una matrice quadrata (di ordine k) che si ottiene scegliendo k righe e k 
colonne di A e considerando gli elementi disposti su tali righe e colonne , 
mantenuti nello stesso ordine. 


Esempio 3.7.2 Data A = 


/ 3 
3 

-2 

\4 


1 0 
2 -1 

1 


1 

3 


1 

2 

6 

0 


1 2 -2 


1\ 
0 
0 
2 / 


se scelgo la l a ,23 e 3 0, 


riga e l a , 3 a e 5 a colonna , ottengo il seguente minore di ordine 3: 


M = 




Chiaramente da ogni matrice si possono estrarre molti minori di vari ordini. 
Se una matrice A è di tipo n x m, il massimo ordine che un minore estraibile 
da A può avere è il minimo tra n e m. Se A è quadrata di ordine n, essa 
stessa è (runico) minore di A di ordine n. 

Definizione 3.7.3 La matrice A si dice avere caratteristica o rango k se: 

1 ) Esiste un minore di A di ordine k con determinante non nullo. 

2) Tutti i minori di A di ordine maggiore di k hanno determinante nullo. 

In tal caso scriviamo car A = k, o anche rk A = k. 

La matrice nulla ha caratteristica zero. 


Esempio 3.7.4 1) Sia 


A = 



(3.8) 
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Osserviamo che il minore 


, ottenuto scegliendo le prime due righe 


2 -1 
3 1 

e le prime due colonne , ha determinante non nullo, mentre tutti i minori di 
ordine 3: 







hanno determinante zero. 
2) Sia 


Quindi la caratteristica di A è 2 


B = 



Poiché det B — — 1^0, car B = 3. 


Si noti che, se A è quadrata di ordine n, allora car A = n se e solo se det ^4^0. 

Si badi bene al fatto che la caratteristica di una matrice è sempre un 
numero intero non negativo, a differenza del determinante che è un numero 
reale. 

Osservazione 3.7.5 Se ho dimostrato che tutti i minori di ordine k hanno 
determinante nullo, è inutile proseguire alla ricerca di minori di ordine an- 
cora maggiore aventi determinante diverso da zero. Infatti, per lo sviluppo 
di Laplace, il determinante di un minore di ordine h si scrive come combina- 
zione lineare di un certo numero di determinanti di ordine k < h e, se questi 
sono tutti nulli, anche il minore di ordine h avrà determinante uguale a zero. 


Osservazione 3.7.6 Le trasformazioni elementari di riga non alterano il 
rango di una matrice, cioè applicando trasformazioni elementari di riga ad 
una matrice si ottengono matrici con lo stesso rango. Quindi matrici equi- 
valenti per riga hanno stesso rango. Lo stesso vale per le colonne. 

Proposizione 3.7.7 Sia A una matrice di caratteristica p e indichiamo con 
Ak xì Ak 2 , . . . , Ak p le colonne che vanno a formare un minore M di ordine 
p con determinante non nullo. Allora ogni altra colonna A^ di A si può 
ottenere come combinazione lineare di queste. Analogamente per le righe. 
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Dim. Senza ledere la generalità, possiamo supporre che le colonne che 
formano M siano le prime p colonne e le righe che formano M siano le prime 
p righe. Sia A' h la colonna che si ottiene considerando i primi p elementi di 
Ah- Sia A la colonna di p elementi così definita: 

A = M~ 1 A ' h . 


Si noti che MA = A' h . 

Se Ài, . . . , X p sono gli elementi di A, sia A'h la colonna così ottenuta: 

A'h = Ah — (Ai Afa + • • • + \ p Ak p ) . 

È facile vedere che i primi p elementi di A'h sono tutti zero. Se A'h è una 
colonna fatta da tutti zeri, allora Ah è combinazione lineare degli A&., come 
volevasi. Se, per assurdo, vi fosse un elemento di A'h diverso da zero, ad 
esempio il j -esimo, consideriamo il minore M' di ordine p + 1 di A ottenuto 
considerando le prime p righe di A più la j-esima e le prime p colonne di A 
più la ù-esima, cioè Ah- Calcoliamo il determinante di M' . Sappiamo che 
il determinante di una matrice è uguale al determinante della matrice che si 
ottiene sottraendo ad una colonna una combinazione lineare delle altre. Nel 
nostro caso, se sottraiamo airultima colonna di M' la combinazione lineare 
delle altre con coefficienti i À^ di cui sopra, si ottiene una matrice M" che ha 
come ultima colonna tutti 0 tranne airultimo posto dove vi è un elemento 
a / 0. Se calcoliamo det M" = det M' secondo l’ultima colonna si ha che: 
detM" = adetM ^ 0. Ma ciò è impossibile, perché allora A avrebbe un 
minore di ordine p + 1, cioè M', con determinante non nullo, mentre, per 
ipotesi, carA = p. 

L’enunciato per le righe si può ottenere ripetendo lo stesso ragionamento 
con t A. □ 


Corollario 3.7.8 Se A' è una matrice ottenuta da A aggiungendo una co- 
lonna B , allora car A < car A' . Inoltre : 

B è una cominazione lineare delle colonne di A se e solo se car A = car A'. 


Determinare la caratteristica di una matrice è un procedimento abba- 
stanza laborioso, perché occorre calcolare molti determinanti. Per ridurre il 
numero di casi da analizzare, si può utilizzare il risultato che segue. 
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Definizione 3.7.9 Se M è un minore di A, un orlato di M è un minore di 
A che ha M come suo minore. 

In sostanza un orlato di M è un minore di A che si ottiene aggiungendo 
a M delle righe e delle colonne di A. 

Proposizione 3.7.10 Sia A una matrice e M un suo minore con det M/ 0 
di ordine h. Se ogni orlato di M di ordine h + 1 ha determinante nullo , allora 
car A = h. 

Dim. Supponiamo che ogni orlato M' di M di ordine h + 1 sia non in- 
vertibile. Basta dimostrare che il rango di A è h. Consideriamo la matrice 
A formata dalle colonne che individuano il minore M. Questa ha rango h. 
Se aggiungiamo a questa matrice una qualsiasi altra colonna Aj di A, per 
l’ipotesi fatta, otteniamo una matrice ancora di rango h. Quindi, per il Co- 
rollarlo 3.7.8, la colonna Aj è combinazione lineare delle colonne di A e, con 
operazione elementari di colonna, può essere rimpiazzata con una colonna 
nulla senza per questo alterare il valore del rango della matrice. Ripetendo 

10 stesso procedimento per ogni colonna di A che non sia tra quelle che in- 
dividuano il minore M, troviamo dia la matrice A ha stesso rango di una 
matrice che come colonne (a meno dell’ ordine) quelle di A oltre a colonne 
nulle, che ha chiaramente rango h. Abbiamo così provato che rk A = h. □ 

Abbiamo dunque il seguente: 

Criterio 3.7.11 (del minore orlato) Una matrice A ha caratteristica r se 
e solo se esiste un minore di A invertibile tale che ogni suo orlato di ordine 
r + 1 non sia invertibile. 

11 risultato sopra permette di semplificare il calcolo del rango di una matrice. 

Ad esempio, nel caso della matrice in (3.8), dopo aver individuato un 

, bastava 

limitarsi a considerare solo i minori di ordine 3 che contenevano M : 

/2 — 1 3 \ /2 -1 1\ 

3 1 2 , 3 14. 

\0 2 - 2 / \0 2 2 / 

Poiché essi hanno determinante nullo, potevamo subito concludere, per il 
criterio dei minori orlati, che car A = 2. 


minore di ordine 2 con determinante non nullo, M = 
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3.8 Esercizi 

1) Provare che, se A e B sono matrici antisimmetriche, allora AB — BA è 
ancora antisimmetrica. 


Determinanti 


2) Calcolare il determinante delle seguenti matrici: 

fi a a 2 \ 

1 b b 2 ; ( 3 . 9 ) 

\l c c 2 ) 


f 0 a b\ 

I —a 0 c ; 
\-b -c 0/ 



fi 

0 

0 a\ 



0 

1 

0 b 

• 


0 

0 

1 c 



\a 

b 

c 0 / 



(k 

1 

0 0\ 


0 

k 

1 0 

• 


0 

0 

k 1 



K 1 

0 

0 k) 


/ 1 


2 

-1 - 

-3\ 

3 


-1 

2 

4 

-2 

-3 

5 

1 

V 4 


1 

-3 

2 / 

(2 

1 

-1 

0 

3 \ 

1 

1 

2 

1 

3 

0 

0 

3 

2 

5 

0 

0 

3 

1 - 

-1 

\0 

0 

0 

0 - 

-v 


( 3 . 10 ) 


( 3 . 11 ) 


( 3 . 12 ) 


( 3 . 13 ) 


( 3 . 14 ) 
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3) Una matrice quadrata della forma 


A = 


(Ai 

0 

• 0 \ 

0 

A2 

... 0 

• • • 

V o 

• • • 

0 

• • • • • • 

• • • aJ 


con Ai , . . . ,A n matrici quadrate (dove 0 indica una matrice nulla di ordine 
opportuno), si dice matrice diagonale a blocchi. Si provi che 


det A = det A\ • • • det A n . 


4) Dimostrare che una matrice reale antisimmetrica di ordine dispari ha 
determinante nullo. Si può dire qualcosa di analogo per quelle di ordine 
pari? 


5) Verificare che se / : 
volte, allora si ha 


R e g : R — > R sono due funzioni derivabili due 


det 


/ r\ 

9 9 


= det 


/ / 
9 9 


H 


H 


6) Provare per induzione il seguente risultato (det. di Vandermonde) 


det 

(1 

1 

• 

• 

• 

Xi 

x 2 

• 

• 

■ 

X 1 “ * 
x 2 • • • 

• • 

• • 

• • 

x n r l \ 

xX 1 

• 

• 

• 

n ì<j{_Xj xi) . 

(3.15) 



X n 

~>2 

X n * * * 

rpLl 1 1 

/ 




Invertibilità di matrici 


7 ) Stabilire se le seguenti matrici sono invertibili e, in caso affermativo, 
calcolare la matrice inversa (verificandone il risultato): 




(3.16) 


(3.17) 
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(2 -1 
3 2 
\0 -2 

/ 2 3 

r 1 2 

\0 1 



•) 



1 


/ 1 0 2 3 \ 

1 - 20-1 
0 3 2 3 

V-i ili/ 

8) Stabilire se esiste una matrice X G M 3 (M) tale che 


1 7 1\ /I 1 1\ /IO 1\ 

1 7 2 X 1 2 3 = 1 1 0 

173/ \1 4 9/ \0 1 1/ 


(3.18) 


(3.19) 


(3.20) 


9) Determinare per quali valori del parametro a G M esiste X G M 3 (M) tale 
che 

/ 2 0 a\ (l 0 0\ /IO -1\ 

-1 1 0X1 -1 0 = a 1 7 . 

\ 1 12/ \1 0 a) \5 0 3 / 


Operazioni di riga - Rango di matrici 


IO) Si riducano con operazioni elementari di riga a forma triangolare le 
seguenti matrici e, nel caso di matrici invertibili, si calcoli l’inversa. 


fi 0 0 

0 1 3 

11-1 

\0 2 0 




0 2 
-2 0 
3 2 

1 1 


3 \ 
-1 
3 
!/ 


(3.21) 


(3.22) 
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11) Calcolare il rango delle seguenti matrici: 


/0 2 — 4 \ 

14-5 
3 1 7 

0 -5 10 
\2 3 0 / 



-2 

9 

7\ 

co 

0 

3 

2 

1 

3 

-1 

2 

V 


Il 1 1 1 7 + A 

12 0 1 0 

4 5 0 4 0 

\1 2 0 1 0 / 


1—3 -1 -2 -1 — 4\ 
2i 0 0 0 0 

2 1110 
12 12 2 

\ 1 2-12 3 / 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


12) Calcolare, al variare del parametro 








Il 0 k 2\ 

11 1 1 ; (3.27) 

\2 k -3 0/ 

Ik - 1 2 \ 

1 -2 0 . (3.28) 

\1 k -2/ 

13) Dimostrare che una matrice reale antisimmetrica ha rango pari. In par- 
ticolare il determinante di una matrice antisimmetrica di ordine dispari è 
nullo. 


14) Data una matrice A di tipo m x n e un suo minore M di ordine p, si 
calcolino quanti sono i minori di A ordine k e gli orlati di M di ordine k. 


Capitolo 4 

SISTEMI LINEARI 


4.1 Sistemi lineari 


Definizione 4.1.1 Si dice sistema lineare di m equazioni in n incognite un 
insieme di m equazioni algebriche tali che siano tutte dei polinomi di pri- 
mo grado in n indeterminate a coefficienti reali (eventualmente complessi) 
uguagliati a zero. 

Un sistema lineare si può scrivere: 


( 

ai\X\ + 0 - 12^2 + * * * + CLl n%n 




a 21 x l H" a 22 x 2 + * * * + &2 n %n 




Uml%l + O j m2%2 + * * * + ttmn^n 

K 



noti. 


La matrice 


A = 


gli aij 

sono 

i coefficienti , e ibi sono i termini 

( a n 

a 12 

. . . O'in ^ 


&21 

• 

&22 

• 

• • • ^2 n 

• 


• 

• 

\^ral 

• 

• 

&m2 

• • 

• • 

. . . a mn J 



che ha al posto (i,j) il coefficiente della j-esima incognita nella i-esima 
equazione, si dice matrice dei coefficienti (o incompleta) del sistema. Essa è 
di tipo m x n, dove m è il numero di equazioni e n è il numero delle incognite. 
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La matrice 



^ a ll 

a 12 

. . . a\ n 

bi\ 

A' = 

a 21 

• 

• 

a 22 

• 

• 

. . . a2n 

• • 

• • 

• 

^2 

• 

• 


• 

^ml 

• 

a m 2 

• 

. . . a mn 

• 

bm J 


che si ottiene da A con raggiunta della colonna dei termini noti, si dice 
matrice completa del sistema. Essa è di tipo m x (n + 1). 

(x A (bi\ 


Se X 


è la n-upla delle incognite e B 


la m-upla dei 


\X n J 


\bnj 


termini noti, il sistema (4.1) si può scrivere in forma matriciale : 


A-X = B. 



Qui A • X è il prodotto riga per colonna tra A , matrice di tipo m x n e X 
matrice di tipo n x 1 il cui risultato è una matrice m x 1 che deve essere 
uguagliata a B, che è dello stesso tipo. 

Definizione 4.1.2 Una n-upla (si, • • • , s n ) di numeri reali è detta una so- 
luzione del sistema dato se, sostituendo alla indeterminata Xi il numero reale 
Si in ogni equazione del sistema, tutte le equazioni risultano verificate. 


Osservazione 4.1.3 Se indichiamo con Ai , . . . , A n le colonne della matrice 

! xi\ 


A, dato X — 


, osserviamo che 


\X n J 


AX — X\A\ + • • • + x n A n , 


cioè il prodotto righe per colonne AX è uguale alla combinazione lineare della 
colonne di A con coefficienti x\, ... ,x n . 

( x i\ 


Dunque, X = 


è soluzione del sistema AX — B se e solo se 



x\Ai + • • • + x n A n — B , 


cioè se e solo se la colonna B dei termini noti è combinazione lineare della 
colonne di A con coefficienti X \, . . . , x n . 
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Il problema consiste ora, dato un sistema lineare, nello stabilire se esistono 
soluzioni e, in caso affermativo, determinarle tutte. 

Nel caso in cui si abbia un sistema di n equazioni in n incognite, per il 
quale quindi la matrice dei coefficienti è quadrata, vale il seguente importante 
risultato: 

Teorema 4.1.4 (Cramer) Sia S un sistema di n equazioni in n incognite 
e sia A la matrice dei coefficienti di S. Se det A à 0, allora il sistema è 
risolubile, ammette una unica soluzione. 

Tale soluzione S è data dalla seguente formula, detta formula di Cramer: 

det Ài 
det A 

dove Ai è la matrice che si ottiene sostituendo in A la i-esima colonna con 
la colonna dei termini noti. 



S (^i, • • • , s n ) •, con Si 


Dim. Poiché det A ^ 0, A è invertibile. Poniamo: 


fh\ 


S = A 1 • B , dove B = 


è la colonna dei termini noti 


Vv 

Allora è chiaro che, scrivendo il sistema in forma matriciale A • X — B , la 
n-upla S sia una soluzione; infatti, sostituita alla n-upla X , si ha che: 


A - S = A - A -1 • B = B . 


Tale soluzione è unica. Siano Si ed S 2 due soluzioni del sistema, sia cioè 
A • Si = B per i — 1,2. Allora A • Si — A • S 2 = B — B = 0, da cui 
A- (Si— S 2 ) = 0. Moltiplicando a sinistra per A~ l si ottiene A -1 -A- (S 1 —S 2 ) = 
(Si — S 2 ) = A -1 • 0 = 0 e quindi Si = S 2 . 

Giustifichiamo ora la formula (4.3). Basta verificare che: 

( 

A~ l B = : 

\s n ) 


con Si = 


det Ài 
det A 


1, . . . , n ; 


quindi la Aesima riga di A 1 moltiplicata per la colonna B dei termini noti 
deve essere uguale a $*. Proviamolo, a titolo di esempio, per i — 1. Per 
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quanto visto nella dimostrazione della Proposizione 3.6.1, la prima riga di 
A _1 moltiplicata per B dà 

b] det Au — b 2 det A 2 1 + • • • + (— 1 ) k+1 b k det A kì + • • • + (— l) 1 +n 6 n det A nì 

det A ’ 

ma il numeratore è proprio det Ài sviluppato secondo la prima colonna. C 


Esempio 4.1.5 Sia dato il sistema S : 


{ 2x + y — z = 1 , 
x + y + 2z = 0 , 
y -\z = 0 . 


2 1-1 

La matrice dei coefficienti è A= \ 1 1 2 ) e si ha det A = — ^ 7 ^ 0. 

0 1 -ì 


Siamo quindi nelle condizioni di applicare il teorema di Cramer. 



det Ai 

det A 


dove Ai 


V = À 2 


2 = Sf dove À 3 







_ 5 _ 
11 ’ 



11 ’ 



2_ 

11 ‘ 


La terna ( jx? — — ji) è soluzione del sistema S . Verificarlo. Tale soluzione 

è unica. 


Definizione 4.1.6 Un sistema lineare si dice omogeneo se tutti i termini 
noti sono zero. 

Osservazione 4.1.7 Ovviamente un sistema omogeneo ammette sempre al- 
meno una soluzione: la n-upla di zeri. 

Se un sistema omogeneo verifica le ipotesi del teorema di Cramer, tale 
soluzione, detta soluzione banale, è anche unica. 
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Purtroppo, il teorema di Cramer non ci è d’aiuto quando: 

1. det A = 0, 

2. A è rettangolare, ossia il numero delle equazioni del sistema è diverso 
dal numero delle incognite, 

essendo A è sempre la matrice dei coefficienti. 

In tutti i casi, però, possiamo applicare il seguente fondamentale risultato: 


Teorema 4.1.8 (Rouché-Capelli) Dato un sistema lineare S, esso è riso- 
lubile se e solo se 

cere A = carA' 

dove A e A' sono rispettivamente la matrice dei coefficienti e la matrice 
completa di S. 

Dim. Per l’Osservazione 4.1.3, il sistema AX = B ha soluzione se e solo se la 
colonna B dei termini noti è esprimibile combinazione lineare delle colonne 
della matrice A. Per il Corollario 3.7.8, questo è possibile se e solo se la 
caratteristica della matrice A e della matrice A' (ottenuta da A aggiungendo 
la colonna B ) coincidono. □ 


4.2 Risoluzione dei sistemi lineari 

A questo punto sappiamo dire se un dato sistema ha soluzioni oppure no. Vo- 
gliamo dare un procedimento algoritmico che permetta di determinare, quan- 
do esistono, tutte le soluzioni del sistema dato, qualora non fosse applicabile 
il teorema di Cramer. 

Schema 1 Consideriamo il sistema in forma matriciale: 

AX = B. 

1) Si verifichi per prima cosa che il sistema abbia soluzioni , cioè per Rouché- 
Capelli , che carA = carA', dove A' è la matrice completa. 

2) Se M è un minore di A (e anche di A') con determinante non nullo di 
ordine massimo, possiamo trascurare, letteralmente cancellandole, tutte le 
equazioni del sistema le cui righe corrispondenti non fanno parte di M . 
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Il sistema così ottenuto è equivalente al sistema dato, cioè i due sistemi 
hanno lo stesso insieme di soluzioni. Infatti se una equazione f(X ) = 0 
si ottiene come combinazione di altre equazioni f(X) = 0 (nel senso che 

f(X) = Xifi(X) H b A nfn(X)), seY è una soluzione comune alVinsieme 

di equazioni fi(X) = 0, ... , fn(X) = 0, Y è soluzione anche dell' insieme di 
funzioni fi(X) = 0, ... , f n (X) = 0, f{X) = 0 e, ovviamente, viceversa. 

3) I termini delle incognite i cui corrispondenti coefficienti formano colonne 
che non stanno in M vengono portati a termine noto e considerati come 
parametri. 

4) Se k è il rango di A, a questo punto abbiamo un sistema di k equazioni 
(abbiamo cancellato tutte le altre) in k incognite, mentre tutte le altre n — k 
incognite sono da considerarsi parametri. Di questo sistema M è la matrice 
dei coefficienti ed ha determinante non nullo. 

5) Si risolve questo sistema con il metodo di Cramer, tenendo conto del 
fatto che i termini noti non sono costanti ma dipendono da n — k parametri, 
parametri che compariranno nella espressione delle soluzioni. 

Esempio 4.2.1 Sia S il seguente sistema: 

{ 2x - \y + w = 1 , 
x+^y + z + w = 1 , 
x — y + w = 1 . 

È un sistema di tre equazioni in quattro incognite. La matrice dei coefficienti 

A è 

(2 0 1\ 

1 | 11- 
Vi -1 0 l) 

Calcoliamo il suo rango. Osserviamo che il minore 



\l -1 0 / 

ha determinante |. Allora car A = 3. 

Necessariamente, essendo una matrice 3x5 ; anche A' ha rango 3, dunque 
il sistema è risolubile. 

Poiché nessuna riga è esclusa dalla formazione di M , non dobbiamo 
cancellare nessuna equazione. 
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L’incognita w, i cui coefficienti formano la quarta colonna di A , esclusa 
da M , è da considerare come parametro e i termini che la contengono vanno 
portati a termine noto. 

Si perviene ad un sistema S' equivalente ad S che ha M coma matrice 

fi — w 

dei coefficienti e come colonna dei termini noti B' = I 1 — w 

\l-w 

Applicando il metodo di risoluzione di Cramer a S' si ha: 


* = St . Mi 

V = St . dove Ah 

* = St . dove M 3 


fi —w 0 \ 

l-u, 1 1 

\l-w -1 0/ 

/ 2 1 -w 0 \ 

1 1-w 1 , 

\1 1-w 0 / 

/2 l-^\ 

1 | 1-w] 

\l -1 1-w) 


da cui x — |(1 — w ) , 


da cui y = — 1(1 — w ) , 


da cui z = (l — w) . 


Al limite, per enfatizzare la trasformazione di w in parametro, cioè la possi- 
bilità di attribuirgli un valore arbitrario e da questo ricavare i corrispondenti 
valori delle altre incognite che vanno a costituire una soluzione, che, ricordia- 
mo, è una n-upla (in questo caso una quaterna) di valori, possiamo introdurre 
come parametro reale t, porre w — t e scrivere : 

Il sistema S ammette come soluzioni tutte le quaterne di numeri reali: 



( x \ 

( U 1 -*) \ 

(\\ 

(A\ 

y 


-la -t) 


2 

3 

+t 

2 

3 

Z 


i ~t 


ì 


-ì 

W 


\ t ) 


\D 


W 


con t Gl. 


Il sistema ammette dunque infinite soluzioni che sono in corrispondenza 
biunivoca con V insieme IR. 

Il numero di parametri che compaiono nella descrizione dell’ insieme delle 
soluzioni di un sistema è n — k dove n è il numero di incognite e k è il rango 

di A. 
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4.3 II metodo di Gauss-Jordan 

Torniamo al caso di un sistema di n equazioni in n incognite in cui la matrice 
dei coefficienti abbia determinente non nullo. In questo caso si può applicare 
il teorema di Cramer ed essere sicuri dell’esistenza di una unica soluzione. 
Inoltre, con la formula data siamo in grado, teoricamente, di determinarla. 
Il problema non teorico ma di calcolo nasce dal fatto che per trovare con 
quel metodo la soluzione occorre calcolare un certo numero di determinanti. 
Fintanto che le dimensioni del sistema sono piccole, va ancora bene, ma non 
appena il sistema diventa molto grande (e in molte applicazioni numeriche 
questa è la regola, avendosi sistemi con migliaia di incognite), il calcolo dei 
determinanti diventa oneroso anche per macchine molto potenti. Vogliamo 
descrivere un metodo, detto di Gauss-Jordan, che permetta di risparmiare 
una considerevole quantità di calcoli. 

Tutto parte dalla seguente semplice osservazione. 

Osservazione 4.3.1 Supponiamo di eseguire le trasformazioni elementari 
di riga sulla matrice completa A' di un sistema. La matrice che si ottiene 
è la matrice completa di un sistema equivalente a quello dato. Infatti , se 
scambiamo due righe di A' , a livello di sistema, significa scambiare V ordine 
delle equazioni. Ma V insieme delle soluzioni, che è V intersezione degli insiemi 
delle soluzioni delle singole equazioni, non cambia (proprietà commutativa 
dell' intersezione) . Per ciò che concerne le trasformazioni di secondo tipo, 
date n equazioni fi(X) = 0, . . . , f n (X) = 0 , Y è soluzione comune a queste 
se e solo se è soluzione comune a fi(X) = 0, . . . , / n _ i(V) = 0 e À 1 / 1 (V) + 
* * * + A n _i/ n _i(X) + fn(X) = 0. 


Schema 2 Dato un sistema di n equazioni in n incognite, sia A la matrice 
dei coefficienti e B la colonna dai termini noti. Supponiamo che det d / 0 
e quindi ci troviamo nelle ipotesi del teorema di Cramer. 

1 ) Ridurre A a forma triangolare mediante trasformazioni elementari di riga 
(assolutamente non di colonna) secondo il procedimento descritto nella di- 
mostrazione di 3.4-8, ottenendo così una matrice A' . 

2) Eseguire contemporaneamente le stesse trasformazioni elementari sulla co- 
lonna dei termini noti, ottenendo una colonna B' . 

3) Poiché la matrice A' ottenuta sopra è triangolare, l'ultima equazione del 
sistema associato A' • X = B' , è di primo grado con una sola incognita e 
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dunque immediatamente risolubile. Quindi : 


x 


K 


n 


a 


nn 


4) Sostituire il valore di x n così ottenuto nella penultima equazione che diven- 
ta allora una equazione nella sola incognita x n -\, che può dunque facilmente 
essere ricavata. 

5) In generale , poiché la matrice dei coefficienti è triangolare , la i-esima 
equazione contiene solo le incognite Xj con j > i. 

Quindi dopo aver ricavato dalle altre equazioni i valori di Xi+ \,...,x n , 
sostituendo questi valori nella i-esima equazione, Punica incognita che resta 
è la X{, che a sua volta viene ricavata. 

Esempio 4.3.2 Sia S il sistema A • X = B , dove : 


A = 


( 0 

i 

0 

1\ 

(l\ 

2 

i 

2 

0 

, B = 

1 

-è -i 

3 

1 

2 

U i 

-2 

V 


w 


1) Poiché an = 0 ? scambio la prima riga di A con Vultima, che ha primo 
termine non nullo. Stessa operazione compio su B. 

Ottengo come nuove matrici: 


/ 1 

ì 

-2 

A 

fl\ 

2 

ì 

2 

0 


1 

ì 

2 

-ì 

3 

1 

5 

2 

\o 

ì 

0 

V 

vJ 


2) Sommo alla terza riga la prima moltiplicata per — 
per B. 

(\ 1 -2 


031 

OH 


1 

2 


Stessa cosa 


2 12 0 


1 

0 2 2 

? 

5 

2 


\0 


1 


3) Sommo alla seconda riga la prima moltiplicata per — fff = —2. Stessa 
cosa per B. 


fi ì 

0 -1 
o 4 

-2 

6 

2 

2 ^ 
-4 

2 

5 

fi\ 

-i 

5 

2 

\0 1 

0 

!/ 

\G 
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A questo punto la prima colonna ha tutti gli elementi nulli tranne il primo. 
Inizio la trasformazione della seconda colonna, come se stessi lavorando sul 

-1 6 -4 

2 


minore 


— i 2 
2 z 

1 0 1 


4) Sommo alla quarta riga la seconda moltiplicata per — ^ = 1. Stessa cosa 
per B. 


(\ 1 -2 2 \ 


fl\ 

0-1 6 -4 


-1 

0 -è 2 2 

5 

5 

2 


\0 0 6 -3/ 


V o / 


5) Sommo alla terza riga la seconda moltiplicata per — 
per B. 

/ 1 1 —2 2 \ / 1 \ 


0-32 

O 22 


— Stessa cosa 


0-1 6 -4 


-1 

0 0-14 

? 

3 

\0 0 6 - 3 ) 


W 


Anche la seconda colonna ha tutti gli elementi situati sotto la diagonale prin- 
cipale uguali a zero. Affrontiamo la trasformazione della terza colonna. 

6) Sommo alla quarta riga la terza moltiplicata per — ^ = 6. Stessa cosa 


per B. 


/I 1 -2 2 


0-1 6 -4 


-1 

0 0-14 


3 

\0 0 0 21 j 


\18/ 


A questo punto la triangolarizzazione di A è completata. 

Iniziamo a ricavare, partendo da x 4 , i valori da attribuire alle incognite 
per avere la soluzione. 

7) X A = 


18 


21 


6 

7* 


3 

7* 


8) —xs + 4-^ = 3 ? da cui x% = 

9) —X 2 + 6 • | — 4 • | = — 1, da cui X 2 

10) X\-\- \ — 2 • | + 2 ■ | = 1, da cui X\ — 0. 
Dunque la soluzione del sistema è 


1 

7 
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4.4 Esercizi 


Sistemi lineari (senza parametri) 

1) Discutere (e quando possibile risolvere) i seguenti lineari 

x + y — 3z = 0 , 

2x-y = l, 

2y - z = 0 ; 

x + y — z = 5 , 

4x — y — 5z = 6 , 
x — 4y — 2z = —4 ; 

2xi - x 2 + 3x 3 + £ 4 = 1 , 

Xi + X 2 - x 3 — 0 , 

X\ — X 4 = 2 : 


' x 1 — x 2 + 2x 3 + x 4 = 1 , 
4xi + x 2 + 2 x 3 = 1 , 

X\ + x 2 + x 3 + x 4 = 0 , 

„ 5xi + 4x 3 + 2x 4 = 3 ; 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


3x + 2y + z = 1 , 

5x + 3y + 3z = 2 , 

7x + 4y + 5z = 3 , 

^x + y — z = 0 ; 

"xi + 2x2 — 3 x 3 + 5x 4 = 3 , 
x\ - x 2 + x 3 - 4x 4 = -5 , 
5x 4 + x 2 — 3x 3 — 2 x 4 = —9 
k 2x 4 + x 2 - 2x 3 + x 4 = -2 . 


(4.8) 


(4.9) 
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Sistemi lineari (con parametri) 


2) Discutere al variare del parametro fceRi seguenti sistemi lineari: 


X\ + x 2 — x 3 + kx± = 0 , 
kx 2 + x 3 - x 4 = 1 , 

3xi + 4x 2 = 2 ; 


(4.10) 


( 2k + 1)ìc + (k + 1 )y + 3 kz = k , 

(2 k - l)x + (k — 2 )y + {2k - T)z = k + 1 , 
3/cx + 2% + (4 k — l)z = 1 ; 

' 2x + y + z = 2 , 

3ic — y + 2z = 6 , 

+ 2y + 3z = 2 , 
k 5x + ky — z = 3/c ; 


" x — 2y + fez = 0 , 



2 re + y = 2 , 

^x — y + z = k . 


(4.11) 


(4.12) 


(4.13) 


3) Discutere, al variare di A ,/iE R, la risolubilità dei seguenti sistemi lineari: 

{ y- Xz = l : 

x — yy + z = X , (4-14) 

liy-Xz = n\ 

{ x + y - z = n, 
y + z = A, 
x + Àz = — y . 


(4.15) 


Capitolo 5 

SPAZI VETTORIALI 


5.1 Spazi vettoriali. Definizione ed esempi 

Definizione 5.1.1 Un insieme V si dice spazio vettoriale sul campo K se 
sono verificate le seguenti condizioni: 

a. Esiste una operazione binaria interna da V x V in V , detta somma, 
indicata con il simbolo +, munito della quale V è un gruppo abeliano. 
Ciò vuol dire che : 
a.l. La somma è associativa, cioè 

u + (v + w) = (u + v) + w , \/u, v, w E V . 

a. 2. Esiste un elemento di V , solitamente indicato con Oy, tale che : 

U + 0y = 0y+U = U, \/u £ V . 

0 y viene detto elemento neutro oppure zero di V . Se ciò non dà luogo 
ad ambiguità, Oy può essere scritto semplicemente come 0. 

a. 3. Per ogni elemento di V esiste un unico elemento, che sarà indicato 
con —u, tale che: 

u + (—u) = (—u) + u = 0 . 

L‘ elemento —u viene detto /’opposto di u. 

a. fi La somma è commutativa, cioè: 

u + v = v + u , Vu, v G V . 
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b. Esiste, inoltre, un 7 operazione daWLxV inV, detta prodotto per scalare, 
indicata nel seguito con il simbolo • (simbolo che in molti casi può essere 
omesso in analogia con la moltiplicazione tra numeri reali), avente le 
seguenti proprietà. 

Per ogni à,/ì6K e per ogni u, v E V , si ha: 
b.l. (À • pi) • u = X • (fi • v), 
b.2. (À + p) • u = À • u + pi • u, 
b.3. \ -{u-\~v) = \- u J r\-v. 

b.f. Se 1 è Videntità del prodotto di¥L, si ha 1 • u = u, per ogni u E V . 

Gli elementi di uno spazio vettoriale in generale si dicono vettori. 

Il campo K viene detto anche campo base e i suoi elementi si dicono 
scalari. 

E sempre importante tener presente il campo base K; per enfatizzare questo, 
specialmente se vi è pericolo di confusione, se E è uno spazio vettoriale su 
K, diremo che V è un K-spazio vettoriale. 

Nel seguito, supporremo sempre K = R e, in qualche caso specificato, 

K = C. 

Facciamo alcune precisazioni. Nella precedente definizione vi è stato un 
abuso di notazione; infatti, il simbolo + è servito per indicare sia la somma 
in V che la somma in K. Queste due operazioni, in alcuni esempi, saranno 
strettamente imparentate, ma in generale tra di esse non vi è alcuna corre- 
lazione. Lo stesso dicasi per il simbolo -, che sta ad indicare sia il prodotto 
tra elementi di K che il prodotto per scalare. 

Il segno • verrà sovente omesso, dopo esserci rassicurati sulla nostra capa- 
cità di percepire i diversi significati indicando, per esercizio, nella definizione 
precedente, quando esso indica il prodotto di K e quando, invece, il prodotto 
per scalare. 

Vediamo ora qualche conseguenza immediata della definizione. 

Osservazione 5.1.2 Sia u un qualsiasi elemento dello spazio vettoriale V 
e À E K uno scalare qualsiasi. Indichiamo con 0 V elemento neutro della 
somma del campo K. Allora si ha: 

1) 0 ' U — Oy . 

Infatti, u = l- u=(l + 0)-u = l- u + 0- u = u + 0-u. 

2) (—1) • u = —u. 

Infatti, Oy = 0 • u = (1 + (—1)) • u = 1 • u + (—1) • u = u + (—1) • u; 
quindi (— 1) • u = — u, per V unicità dell 7 opposto di u. 
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3 ) A • Oy — 0 y. 

Infatti, AOy = A(0y + Oy) = AOy -|- AOy. 

Esempio 5.1.3 1) R. 

L’insieme dei numeri reali è uno spazio vettoriale su se stesso, considerando 
come operazione di somma la stessa somma tra numeri reali e come prodotto 
per scalare il prodotto stesso di R. Lasciamo come esercizio la verifica in 
dettaglio. 

2) R 2 . 

Consideriamo come somma tra coppie la seguente operazione: 

(ai, 6i) + (a2, 62) = ( a i + a 2, b\ + ^2) 

e come prodotto per scalare il seguente: 

A (a, b) = (Aa, A6). 

Anche in questo caso è elementare verificare che sono soddisfatte tutte le pro- 
prietà della definizione di spazio vettoriale. Procedendo in modo analogo, si 
può definire una struttura di R -spazio vettoriale su R n . 

3) M m , n (R). 

L’insieme delle matrici dello stello stesso tipo a valori reali forma un R- 
spazio vettoriale con le usuali operazioni di somma tra matrici e di prodotto 
per uno scalare viste a loro tempo. 

4) C. 

L’insieme dei numeri complessi è un C -spazio vettoriale con le stesse opera- 
zioni di campo di cui è munito, in completa analogia con l’esempio 1). Non 
è però ozioso notare che è anche un R -spazio vettoriale. 

5) R[x]. 

L’insieme dei polinomi a coefficienti reali, con le solite operazioni di somma 
e prodotto per un numero reale. 

6 ) JXAR). 

L’insieme i cui elementi sono le funzioni f : A — > R ? dove A è un insie- 
me qualsiasi. Per ogni coppia di elementi f,g di ^(A,R) definiamo allora 
h = / + g ponendo: 

i) h(x) = f(x) + g{x) \/x G R. 

e definiamo il prodotto per scalare (A/) ponendo: 

ii) (A f)(x) = A f(x) \/x G R. 

Teniamo presente che f{x), g{x), A sono numeri reali e quindi le operazioni 
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a secondo membro nelle due uguaglianze sono tra elementi di R. La verifica 
dettagliata del fatto che M), munito delle operazioni suddette, sia uno 
spazio vettoriale reale è lasciata per esercizio. 

7) L’insieme dei vettori geometrici. 

Quello che segue è un esempio importante per visualizzare intuitivamente gli 
spazi vettoriali. Cominciamo richiamando alcune definizioni. Per segmen- 
to orientato si intende un segmento nello spazio ordinario della geometria 
di Euclide (quello studiato nelle scuole superiori, per intendersi) del quale 
si sia stabilito quale sia il primo estremo e quale il secondo. Ovverosia, se i 

punti A e B sono rispettivamente il primo ed il secondo estremo del segmento 

— > 

AB, tale segmento orientato verrà indicato con AB , dove la freccia sta ad 

evidenziare la relazione d’ordine tra gli estremi del segmento. Chiaramen- 

< — 

te il simbolo AB sta ad indicare lo stesso segmento ma con la orientazione 
opposta degli estremi, quindi, in definitiva, un diverso segmento orientato. 
Poniamo ora, sull’insieme di tutti i segmenti orientati la seguente relazione 
di equivalenza : 

— » — > 

Due segmenti orientati AB e CD sono equivalenti se e solo se essi giac- 
ciono sullo stesso piano ed inoltre il quadrilatero ABDC , avente come lati i 
segmenti AB, BD , DC, CA , è un parallelogrammo. 

In altre parole, due segmenti orientati sono equivalenti, secondo la definizio- 
ne, se giacciono su due rette parallele e traslando una di queste rette sempre 
parallela a se stessa, è possibile portare i due segmenti a coincidere in guisa 
tale che il primo ed il secondo estremo di uno vada a sovrapporsi, rispetti- 
vamente, al primo ed al secondo estremo dell’altro. Debbono perciò avere la 

V 

stessa lunghezza, la stessa direzione e lo stesso verso. E per questo facile 
convincersi che la precedente relazione è effettivamente di equivalenza. 

Si dice vettore geometrico una qualsiasi classe di equivalenza della pre- 
cedente relazione. 

Si dice vettore applicato nel punto P, il rappresentante di una classe di 
equivalenza avente come primo estremo il punto P. 

Definiamo sull’insieme V dei vettori geometrici una operazione di somma ed 
una operazione di prodotto di un vettore per un numero reale che daranno a 
V la struttura di spazio vettoriale. 

La somma. Dati due vettori geometrici u e v, cioè due classi di equiva- 

— > 

lenza, sia AB un rappresentante di u. Si scelga allora, come rappresentante 

— > 

di v, il vettore applicato nel punto B . Sia esso il segmento BC. Allora, per 
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definizione, la loro somma u + v sarà la classe di equivalenza del segmento 

AC. 

y 

In primo luogo , se A coincide con B , allora il segmento orientato AA 
rappresenta il vettore nullo e risulterà lo zero di V. 

Semplici costruzioni geometriche permettono di verificare che tale opera- 
zione gode delle proprietà commutativa e associativa. 

y 

Per inciso, se AB è un rappresentante di w, il vettore — v è rappresentato 

y 

dal segmento orientato BA. 

Il prodotto per scalare. Se v è un vettore e k è un numero reale, allora 
k\ sarà il vettore avente la stessa direzione di v, lunghezza uguale a quella 
di v moltiplicata per \k\ e stessa direzione di v, se k > 0, oppure direzione 
opposta, se k < 0; infine, se k = 0, allora kv è il vettore nullo. 

Con semplici ma lunghe costruzioni geometriche è possibile verificare che 
tutte le proprietà richieste per essere uno spazio vettoriale sono soddisfatte. 


5.2 Sottospazi vettoriali 

Sia V uno spazio vettoriale sul campo K, e sia W un sottoinsieme di V. Dati 
due elementi di W, è possibile sommarli e moltiplicarli per degli scalari, ma 
in generale i risultati di tali operazioni non è detto che appartengano a W. 
Ad esempio, sia C = {(x, y) G R 2 1 x 2 + y 1 = 1}. C è un sottoinsieme di R 2 
che contiene, ad esempio, (1, 0) e (0, 1) ma non contiene la loro somma (1, 1) 
né il prodotto 3 • (1, 0) = (3, 0). 

Definizione 5.2.1 Un sottoinsieme non vuoto W di uno spazio vettoriale 
V è un sottospazio vettoriale di V, se W è chiuso rispetto alle operazioni di 
V, ossia se: 

a) La somma di due elementi di W appartiene ancora a W ; 

b) Il prodotto di qualsiasi elemento di W per uno scalare appartiene ancora 
a W. 


Osserviamo che l’insieme C dell’esempio precedente non è un sottospazio 
vettoriale di R 2 , mentre il sottoinsieme D = {(x,y) G R 2 1 x = y} lo è. 
Infatti dati due elementi di D , i.e. (a, a) e (6, ò), la loro somma (a + b,a + b) 
appartiene ancora a D e se k G R, allora k(a, a) = ( ka , ka) appartiene ancora 
a D. 
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Osservazione 5.2.2 Se W è un sottospazio di V, allora a sua volta W è 
uno spazio vettoriale. La verifica di questo fatto viene lasciata per esercizio. 
In particolare , 0 V G W. 


Osservazione 5.2.3 Le condizioni a) ab) della definizione 5.2.1 sono indi- 
pendenti e debbono essere soddisfatte entrambe affinché W sia effettivamente 
un sottospazio vettoriale. Infatti consideriamo i seguenti esempi: 

a) ZcR è chiuso rispetto alla somma (è un sottogruppo di R/ ; ma non certo 
rispetto al prodotto per scalare . 

b) A — {(x,y) G R 2 1 xy = 0}, il sottoinsieme di R 2 costituito dalla unio- 
ne degli assi cartesiani. A è chiuso rispetto al prodotto per scalare; infatti 
(x,y) G A se e solo se x — 0 oppure y — 0. Allora , se (x,y) G A, anche 
k(x , y) = ( kx , ky) ha almeno una componente nulla e quindi appartiene ad A. 
Ma A non è chiuso rispetto alla somma. Infatti, (1,0) e (0, 1) appartengono 
ad A, ma non la loro somma (1, 1). 


Esempio 5.2.4 Sia S = {X G R n | A-X = 0} ; dove A è una matrice mxn e 
0 è lo zero di R m ; cioè V insieme delle soluzioni del sistema omogeneo avente 

V 

A come matrice dei coefficienti. E facile verificare che S è un sottospazio 
vettoriale di R n . 


5.3 Dipendenza lineare 

Dati alcuni elementi di uno spazio vettoriale V su K, le operazioni di cui V 
è munito ci permettono di costruirne degli altri nella maniera che segue: 

Definizione 5.3.1 Un elemento v G V si dice combinazione lineare degli 
elementi ... ,v n di V se risulta che v = Ài^i + • • • + \ n v n V er opportuni 
coefficienti appartenenti al campo base K. 

Esempio 5.3.2 II vettore (1,0,2) G R 3 è combinazione lineare di (1/2, 1, 1) 
e di (0, 1/3,0) con coefficienti 2 e —6. Infatti: 

(1, 0,2) = 2(1/2, l,l) + (-6)(0, 1/3,0). 

Proposizione 5.3.3 Siano Vi,...,v n elementi di uno spazio vettoriale V. 
Consideriamo il seguente sottoinsieme di V : 


Vq = {v G V | V = Ài^! + * * * + \ n Vn , con \ £ K} , 
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cioè Vq è l’insieme di tutte le combinazioni lineari dei vettori rq, . . . , v n . Allo- 
ra Vq è un sottospazio vettoriale di V . Esso si indica anche con C{v 1? . . . , v n }. 

Dim. Siano u e w due elementi di Vq. Per la definizione di Vq, risulta che 

u = À!Ui + - • • + X n v n ew = HiVi~\ \~ia n v n con A* e /q opportuni coefficienti. 

Allora il vettore somma u + w risulterà: 

U + W = AiUi H h \ n V n + LlVl H h HnV n 

— (Al + /i 1)^1 + * * * + (A n + t 1 n) v n • 

Quindi anche u + w si esprime come combinazione lineare dei vettori Vi e 
dunque appartiene a Vq. 

Risulta così soddisfatta la proprietà a) della definizione 5.2.1 , Procediamo 
a dimostrare che vale anche la b). 

Se u £ Vq e k è uno scalare, e se u = Aitq + • • • + A allora ku = 
kXiVi + • • • + k\ n v n) che quindi è, per definizione, elemento di Vq. □ 

Osserviamo che lo zero di V è sempre esprimibile come combinazione lineare 
di vettori: basta infatti prendere i coefficienti della combinazione tutti nulli. 


Definizione 5.3.4 Dato un sottoinsieme X di V , il sottospazio vettoriale 


A 

X = {v G V | v = A 1 U 1 + • • • + A n v n , con A* e K , Vi e X} 


viene detto sottospazio generato da X . Nulla vieta che coincida con V stesso. 

Un sottoinsieme X di V forma un insieme di generatori di V , se accade 
che : 


V = {ve V 


v 


Ai^i + • • • + A n v n , con 


In altre parole, se lo spazio X coincide con tutto V. 

Esempio 5.3.5 / vettori (1,1), (0,2) e (—1,1) sono dei generatori di R 2 . 
Infatti, sia (a, b) un qualsiasi elemento di R 2 . Voglio provare che esso si 
ottiene come combinazione lineare dei tre vettori dati, cioè che esistono tre 
coefficienti reali x,y,z, tali che : 


(a, b) — x(l, 1) + j/(0, 2) + z(—l, 1) . 


(5.1) 
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Utilizzando le operazioni di spazio vettoriale di R 2 , possiamo scrivere : 

(a, b) = (x, x ) + (0, 2 y) + (—z, z) = (x - z, x + 2y + z ) . 

Questa uguaglianza tra vettori implica che debba sussistere V uguaglianza tra 
le due componenti e perciò risulta essere equivalente al seguente sistema: 


x — z — a , 
x + 2y + z = b . 


La matrice dei coefficienti di tale sistema è 


1 0 
1 2 



che ha rango 2. 


Per il teorema di Rouché- Capelli, il sistema è sempre risolubile. Quindi, 
comunque presi a e b, troverò sempre degli opportuni coefficienti x,y,z tali 
che valga la relazione (5.1). 


Supponiamo che v G V si possa esprimere come combinazione dei vettori 
Vi, . . . ,v n , cioè esistono dei coefficienti À* tali che v = À^i + • • • + \ n v n . 
Questi coefficienti sono unici oppure possono esistere degli altri coefficienti 
Pi diversi dai À* ma tali che ugualmente si abbia v = p\V\ + • • • + p n v n ? 
Questa domanda motiva la seguente definizione. 

Definizione 5.3.6 Dei vettori v\, . . . ,v n si dicono linearmente dipendenti, 
se esistono dei coefficienti Xi, i = 1 , . . . , n, non tutti nulli, tali che 


Ài^i + • • • + \ n v n — 0v • 


In caso contrario, si dicono linearmente indipendenti. 

Un sottoinsieme A C V è detto insieme linearmente dipendente se con- 
tiene un numero finito di elementi ai, ... linearmente dipendenti. In caso 
contrario, se tutte le n-uple di elementi di A sono linearmente indipendenti, 
allora A viene detto insieme linearmente indipendente. 

Esempio 5.3.7 I vettori (1,1), (0.2) e (—1,1) di IR 2 , sono linearmente di- 
pendenti. Infatti supponiamo che sia : x(l, 1) + y( 0,2) + z(— 1, 1) = (0,0). 
Sviluppando i calcoli si ha il seguente sistema: 

f x — z — 0 , 
x + 2y + z = 0 . 
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Questo sistema omogeneo ha, oltre la soluzione nulla, infinite soluzioni non 
banali, ad esempio: 



Quindi la combinazione lineare con tali coefficienti non nulli è ancora (0,0). 


so- 


Esempio 5.3.8 Le matrici A = ( J ^ , B = ^ , C = ^ ^ J 

no linearmente indipendenti. Infatti, siano X , ji ev i coeff denti di una com- 
binazione lineare di A, B e C che dà la matrice nulla ( ^ ^ . Sviluppando 


0 0 


i conti, si ottiene il seguente sistema: 



Le soluzioni di questo sistema sono i coeff denti di tutte le possibili combi- 
nazioni lineari delle tre matrici date che danno la matrice nulla. La matrice 
dei coeff denti di questo sistema è: 


M = 



det(M) = 2, quindi per il teorema di Cramer la sua unica soluzione è quella 
banale. Si può perciò concludere che Vunica combinazione lineare di A, B , 
C che dà la matrice nulla è quella con i coeff denti tutti nulli. Le tre matrici 
sono dunque linearmente indipendenti. 


Osservazione 5.3.9 La nozione di vettori linearmente indipendenti ci per- 
mette di rispondere al quesito posto più sopra. Supponiamo che il vetto- 
re v sia combinazione lineare dei vettori vi,...,v n con coeff denti Xi, cioè 
v = À^! + • • • + X n v n . Supponiamo anche che esistano degli altri coeff denti 
Ili tali che ugualmente si abbia v = pivi + • • • + ti n v n . Allora abbiamo che : 


(Ai^i H b A n v n ) - (/ìi^i H b fi n v n ) = v - v = 0 V . 
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Raccogliendo i termini del primo membro, si ha: 

(Ài P>l)Vl * * * ~\~ (A n t^n)^n Oy . 

Questa è una combinazione lineare dei vettori Vi che dà lo zero di V. Dunque, 
se i vettori V{ sono linearmente indipendenti, tutti i coefficienti della combi- 
nazione debbono essere nulli; ciò implica che, per ciascun indice i, Xi = /i* 
e quindi vi è la unicità dei coefficienti della combinazione dei vettori Vi che 
esprime v. 

Viceversa, se i vettori Vi fossero linearmente dipendenti, per definizione, 
esisterebbero dei coefficienti non tutti nulli tali che si abbia: u\Vi + • • • + 
VnVn — Oy. Se v — Ài^i + • • • + X n v n , allora si ha: 

V = V-\-§y = + * • •+\ n V n + l'iVi + - * -+U n V n = (Xi + 1Si)Vi + - • • + (X n +v n )v n . 

Quindi anche la combinazione lineare dei vettori vi,...,v n con coefficienti 
hi = Xi + Vi esprime il vettore v; ma, poiché i Vi non sono tutti nulli, qual- 
cuno dei ki sarà diverso dal corrispondente Ày In definitiva, due diverse 
combinazioni lineari dei vettori Vi danno lo stesso vettore v. 


Osservazione 5.3.10 Un sottoinsieme di un insieme di vettori linearmente 
indipendenti è linearmente indipendente. 

Proposizione 5.3.11 Dei vettori vi , . . . , v n , sono linearmente dipendenti se 
e solo se uno di essi può scriversi come combinazione lineare degli altri. 


Dim. Poiché i vettori dati sono linearmente dipendenti, esiste una loro 
combinazione lineare con coefficienti non tutti nulli che dà lo zero di V. Sia 
Ài^i + • • • + X n v n = Oy, con uno dei Xi 0, e.g. Ài 0. 

Allora, dividendo i membri dell’ uguaglianza per —X\ e sommando poi 
si ottiene: 



Viceversa, se v 1 = À 2 r ? 2 H b X n v n , allora v 1 — X 2 v 2 — X n v n = Oy, che è 

una combinazione lineare che dà lo zero e avente i coefficienti non tutti nulli, 
essendo il primo coefficiente 1. □ 


Corollario 5.3.12 Siano . . . ,v n vettori linearmente indipendenti di uno 
spazio vettoriale V. Dato v G V, si ha che . . . ,v n ,v sono linearmente 
indipendenti se e solo se v non è combinazione lineare di v \, ... ,v n . 


5.4. BASI E DIMENSIONE DI UNO SPAZIO VETTORIALE 
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5.4 Basi e dimensione di uno spazio vettoriale 

Definizione 5.4.1 Un sottoinsieme B di uno spazio vettoriale V è una base 
per V , se B è un insieme di generatori di V linearmente indipendenti. 


Esempio 5.4.2 I vettori (1,1,1), (1,— 1,0) e (1,0,0) costituiscono una base 
perR 3 . Verifichiamo che sono linearmente indipendenti. Consideriamo una 
loro combinazione lineare che dia lo zero di R 3 . Cioè : 

2 ( 1 , 0 , 0) = (0,0,0); questa, una volta eseguite le operazioni e uguagliate le 
componenti dei vettori, è equivalente al sistema omogeneo: 

f x + y + z = 0 , 

x-y = 0, 

x = 0 . 


La matrice dei coefficienti di questo sistema è 



Il suo determinante è 1, perciò il sistema ha, come unica soluzione, la solu- 
zione banale; quindi i vettori sono linearmente indipendenti. 

Dimostriamo che essi sono anche dei generatori. Sia (a, 6, c) un generico 
elemento di R 3 . I coefficienti di una combinazione lineare dei vettori (1, 1, 1), 
(1,— 1,1) e (1,0,0) che dia (a, 6, c) sono le soluzioni del seguente sistema: 


x + y + 2 = a , 
x-y = b, 


X — c. 


Per il teorema di Cramer, esso è sempre risolubile e dunque, qualsiasi ele- 
mento di R 3 è esprimibile come combinazione lineare dei tre vettori dati. 


Osservazione 5.4.3 II procedimento precedente può essere generalizzato nel 
modo seguente: 

dati k vettori di R n ; sia M la matrice ottenuta scrivendo sulla i-esima 
colonna le coordinate dello i-esimo vettore. M risulta così essere una matrice 
n x k. I vettori dati sono linearmente indipendenti se e solo se M ha rango 
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esattamente k. Infatti, i coefficienti di una combinazione lineare dei vettori 
dati che dia lo zero sono le soluzioni del sistema omogeneo avente come 
matrice dei coefficienti proprio la matrice M . Per il teorema di Rouché- 
Capelli, Lunica soluzione di tale sistema è quella nulla se e solo se M ha 
rango k. 

Esempio 5.4.4 L'insieme dei monomi {1, x, x 2 , x 3 , . . . , x n , . . .} è una base 
dello spazio dei polinomi con coefficienti reali R[x], Infatti ogni polinomio, 
per definizione, si scrive come combinazione lineare di un sottoinsieme finito 
opportuno di tale insieme. Inoltre è ovvio verificare la indipendenza lineare. 
Una combinazione lineare degli x n è il polinomio nullo se e solo se tutti i 
coefficienti sono nulli. 


Esempio 5.4.5 II numero reale 1 è una base per 
vettoriale. 


inteso come R -spazio 


In generale, uno spazio vettoriale può avere più di una base, come mostra 
l’esempio di R 3 , che ammette come basi, come visto in precedenza, i vettori 
(1,1,1), (1,1,0) e (1,0,0), ma anche, come è facile verificare, l’insieme dei 
vettori 

e i = (1, 0, 0) , e 2 = (0, 1, 0) , e 3 = (0, 0, 1) , 
detta base canonica di R 3 . 

Più complicato è stabilire se uno spazio vettoriale abbia una base ed 
inoltre sarebbe interessante stabilire quale relazione intercorre tre due diverse 
basi di uno stesso spazio vettoriale. Il seguente teorema, di cui omettiamo la 
dimostrazione, dà una risposta a questi interrogativi. 

Teorema 5.4.6 Ogni spazio vettoriale ammette almeno una base e due di- 
verse basi di uno stesso spazio vettoriale sono equipotenti. 

Definizione 5.4.7 Si dice dimensione di uno spazio vettoriale il numero di 
elementi di una sua base. In particolare, uno spazio vettoriale è detto di 
dimensione finita se una sua base ha un numero finito di elementi. 

Lo spazio vettoriale {0} ha dimensione zero. 

La precedente definizione è ben posta in virtù del teorema 5.4.6. Ad esempio, 
R 3 ha dimensione 3 perché abbiamo trovato una base con tre elementi. 
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Osservazione 5.4.8 1) In generale, R n ha dimensione n. Infatti, in analo- 
gia con UEsempio 5.4-2, i vettori = (0, . . . , 0, 1, 0, . . . , 0), che hanno tutte 
le coordinate nulle tranne la i- esima, che invece vale 1, costituiscono una 
base di R n che ha esattamente n elementi. 

2) Lo spazio dei polinomi R[ir], invece, ha dimensione infinita. Infatti la sua 
base indicata neW Esempio 5.4-4 ha infiniti elementi. 

3) Lo spazio M m?n (R) ha dimensione nm. Una sua base è data dalle matrici 
Eij, le matrici che hanno 1 al posto (i,j) e 0 altrove. 

In particolare, lo spazio M n (R) delle matrici quadrate di ordine n ha 
dimensione n 2 . 

Osservazione 5.4.9 La dimensione di uno spazio vettoriale dipende anche 
dal campo base. Ad esempio V insieme dei numeri complessi ha dimensione 1 
se viene considerato come spazio vettoriale complesso, mentre ha dimensio- 
ne 2 come spazio vettoriale sui numeri reali (una base di C come R -spazio 
vettoriale è V insieme {1 ,i}). 

Proposizione 5.4.10 Un sottoinsieme {òi, . . . , b n } di uno spazio vettoriale 
V di dimensione finita, costituisce una base per V se e solo se ogni elemento 
v di V si può scrivere, in modo unico, come combinazione lineare dei vettori 
della base bi. 

Dim. Supponiamo che ò n } sia una base di V. Poiché sono dei 

generatori, per definizione, ogni vettore v di V si può scrivere come loro 
combinazione lineare. Inoltre, poiché sono linearmente indipendenti, i coeffi- 
cienti della combinazione lineare dei bi che esprime ciascun v sono, per quanto 
visto nell’ Osservazione 5.3.9, unici. 

Viceversa, se ogni vettore v può essere scritto come combinazione lineare 
dei &i, . . . , 6 n , essi sono, per definizione, dei generatori di V. Inoltre, siccome 
per ipotesi Punica combinazione dei bi che dà lo zero è quella con i coefficienti 
tutti nulli, essi sono anche linearmente indipendenti. Pertanto {òi , . . . , b n } è 
una base di V . □ 

Teorema 5.4.11 (Completamento ad una base) Sia dato uno spazio V 
di dimensione finita. Se tq, . . . , Vh sono dei vettori di V linearmente indipen- 
denti, allora esistono dei vettori uq, . . . , Wk tali che {v \, . . . , Vh, Wi , . . . , Wk} 
sia una base di V. 

I vettori Wi, ... si possono prendere appartenenti ad un insieme as- 
segnato di generatori di V. 
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Dim. Dimostriamolo per induzione su h. Se h = 0, basta prendere come 
completamento una base qualsiasi di V. Supponiamo che l’enunciato sia vero 
per h; proviamo che vale per h + 1. Siano rq, . . . , v k +i dei vettori indipendenti 
di V e siano allora wi, . . . ,Wk i vettori che completano tq, . . . , v k ad una base 
di V, che esistono per ipotesi induttiva. (Infatti, per l’Osservazione 5.3.10, 
Vi ,...,Vh sono anche essi indipendenti). Allora v^+i si può scrivere come 
combinazione lineare dei rq, . . . , vq, . . . , Wk- Scriviamo: 

v h+ i = XiVi H b A h Vh + HiWi H b /ifcWfc • (5.2) 

Almeno uno dei è diverso da zero, altrimenti Vh + 1 si scriverebbe come 
combinazione dei soli v \ , . . . , rq* e risulterebbero perciò i vettori rq, . . . , v k +i 
linearmente dipendenti. Per comodità, supponiamo che /i*. ^ 0. Sottraendo 
da entrambi i membri Vh+ 1 , dividendo per —fi k ed infine sommando vq., alla 
fine si perviene alla seguente scrittura: 

Wk = kivi H b k h+1 v h+ i + U 1 W 1 H b iSk-iWk-i ■ (5.3) 

Ovverosia, w k si scrive come combinazione di rq, . . . , v^+i e vq, . . . , Wk-i- 
Dimostriamo che {iq, . . . , Vh+u vq, . . . , w k _i} è una base di V e abbiamo 
finito. 

1) {tq, . . . , Vh+i,wi , . . . , Wk- 1 } sono dei generatori di V. 

Infatti, sia v E V] v = t\Vi-\ \-t h v k + s ivq + - • • + s k w k , perché, ricordiamo, 

Vi, . . . , Vh, vq, • • • , Wk sono una base di V. Ma allora, sostituendo in questa 
espressione a w k la combinazione (5.3), si ha: 

v = tivO HhVh+siwO bs/c(&qviH \-k h +iVh+i+viW\-\ \-v k -iW k -i)- 

Quindi ogni v si scrive come combinazione dei . . . , v^+i, rei, ... , w k -\. 

2) {v 1? . . . , Vh+uWi , . . . , uq_i} sono linearmente indipendenti. Supponiamo, 

per assurdo, che: 0 V = a x v i + • • • + + &ivq + • • • + ma i 

coefficienti e bi non siano tutti nulli. In particolare, non sarà nullo a^+i, il 
coefficiente di v^+i, perché, se lo fosse, sarebbero nulli anche tutti gli altri, in 
quanto rq, . . . , v kl vq, . . . , uq_i sono linearmente indipendenti. Sostituendo 
in questa espressione, la combinazione (5.2) qui sopra, si ha: 

0y = aiVi+’ • -+CLh+ i(AiVi+* • •+A/ l V/ l +/i.ivq+- * ’P^ k w k )+biWi+- • -+bk~iWk-i . 

Questa è una combinazione dei rq, . . . , rq, vq, . . . , w k che dà lo zero. Il coef- 
ficiente di w k è dh+i/ik che è non nullo perché entrambi i fattori sono diversi 
da zero. Ciò è assurdo perché tq, . . . , v^, vq, . . . , w k costituiscono una base di 
V e perciò sono linearmente indipendenti. □ 
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Questo teorema ha molte conseguenze, utili da adoperare nelle applicazio- 
ni. Facciamone una lista. Di alcune proprietà daremo una dimostrazione; 
dimostrare le altre viene lasciato come esercizio teorico. 

Teorema 5.4.12 Da un insieme finito di generatori dello spazio vettoriale 
V si può sempre estrarre una base di V , cioè un sottinsieme che sia una base 
diV. 

Teorema 5.4.13 Sia V uno spazio vettoriale V di dimensione n. Dati n 
vettori di V , essi sono lineamente indipendenti se e solo se sono un insieme 
di generatori di V. 

Corollario 5.4.14 Sia V uno spazio vettoriale di dimensione n. SeW è un 
sottospazio vettoriale di V tale che dimV = dim W = n, allora V = W. 

5.5 Coordinate 

Sia B una base dello spazio vettoriale V , di dimensione n. Ordiniamola, di 
modo che sia ben specificato quale sia il primo elemento, il secondo e così 
via. 

Definizione 5.5.1 Sia B = {òi , . . . ,ò n } una base ordinata dello spazio vet- 
toriale V . Si definiscono coordinate del vettore v E V , rispetto alla base B, 
la n-upla di scalari (aq, . . . , x n ) tali che v = X\b\ + • • • + x n b n . 

Facciamo vari commenti. 

In primo luogo questa definizione è ben posta: ad ogni v E V risulta 
assegnata in modo univoco la n-upla delle sue coordinate rispetto ad una 
base. Infatti, in virtù della Proposizione 5.4.10, esiste sempre una combina- 
zione lineare di elementi della base data che esprime v e i coefficienti di tale 
combinazione sono unici. 

In secondo luogo, le coordinate di un vettore dipendono strettamente, 
come è ovvio, dalla base scelta e anche dalhordine con cui vengono considerati 
gli elementi della base. Più sotto vedremo in estremo dettaglio le regole di 
questa dipendenza. Ad ogni modo, qualunque sia la base, lo zero di V 
avrà sempre come coordinate la n-upla (0, ... ,0). Poiché possiamo scrivere 
b\ — 1 • b\ + 0 • 62 + * * * + 0 • b n ed inoltre tale scrittura per esprimere b\ come 
combinazione dei bi è unica, le coordinate di 6 1? nella base { 6 1? . . . , ò n }, sono 
(1,0,. . . ,0). Più in generale, l’i-esimo vettore della base, b^ ha coordinate 


86 


CAPITOLO 5. SPAZI VETTORIALI 


e* = (0, . . . , 0, 1, 0, . . . , 0), la n-upla che ha tutte le componenti nulle tranne 
la i-esima, che vale 1. 

Osservazione 5.5.2 Dalla definizione precedente è immediato dedurre che 
le coordinate della somma di due vettori v e w sono la somma, componente 
per componente, delle coordinate rispettive di v e w e che le coordinate di Xv 
sono le coordinate di v moltiplicate per X. 

Poniamoci dunque la questione fondamentale: supponiamo di conoscere le 
coordinate di un vettore v in una data base B. Quali saranno le coordinate 
di v in un’altra base £>'? Per rispondere, consideriamo alcune notazioni. 

1) B — {&i, . . . , b n } verrà indicata come la vecchia base ; B' = {b [, . . . , b' n } 
verrà indicata come la nuova base ; 

2) Supponiamo di conoscere, per ogni indice j, con j = 1, . . . , n, i coefficienti 
delle combinazioni lineari: bj = aijb[ + • • • + a n fb' n . 

Il termine è il coefficiente di b[ nella combinazione che esprime il j- 
esimo elemento bj della vecchia base. In altri termini, il vettore (ay, . . . , a n j) 
rappresenta le coordinate del j-esimo vettore della vecchia base rispetto alla 
nuova base B' . 

3) Sia M|, la matrice (cqQ che ha, cioè, come vettori colonna le n-uple 
(aij, . . . , a n j) delle coordinate dei vettori della vecchia base rispetto alla nuova 
base. La matrice viene detta matrice di cambiamento di base o anche matrice 
di passaggio dalla base B alla base £>'; essa è una matrice quadrata di ordine 
n. 


Proposizione 5.5.3 Conoscendo la matrice tra due basi B e B' e le 
coordinate (aq,...,a; n ) di un vettore v rispetto alla base B, le coordinate 
(xf . . . , x f n ) di v rispetto alla nuova base si determinano nel modo seguente: 


^ x\ ^ 



• 

0 

0 

= Mg/ 

0 

0 

0 

V'J 


\x n J 



Dim. Dire che le coordinate del vettore v rispetto alla vecchia base B 
sono (aq, . . . , x n ) significa che v = aq b\ + • • • + x n b n . Per come è definita la 
matrice di passaggio si ha anche che bj = aijb[ + . . . + a n jb' n . Effettuando le 
sostituzioni opportune, si ottiene che: 


v 


Xi(anb[ + • • • + a n ib' n ) + • • • + x n (ai n b\ + • • • + a nn b ' n ) - 
(aqan + • • • + x n a\ n )b f -y + • • • + (aqa n i + • • • + x n a nn )b' n . 
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Osserviamo che il coefficiente del vettore della nuova base b[ è il prodotto della 
i-esima riga della matrice di passaggio e la colonna delle vecchie coordinate 
di v. □ 


Esempio 5.5.4 Consideriamo i tre seguenti vettori in R 3 , bi = (1,1,0), 
62 = (0,0, 1) e bs = (—1,0, 1). Essi sono una base di R 3 . Infatti la matrice 
A che ha i tre vettori dati come colonne , 






ha determinante non nullo. Quindi i tre vettori sono linearmente indipen- 
denti e, in virtù del Teorema 5.4-13, costituiscono una base di R 3 . 

Vogliamo determinare la matrice di passaggio dalla base canonica 
C = {(1, 0, 0), (0, 1, 0), (0, 0, 1)} alla base B. Per prima cosa determiniamo le 
nuove coordinate dei vettori della vecchia base. Siano (a 11? a 2 i, a 31 ) le nuove 
coordinate del vettore (1,0,0). Dovrà essere verificata la seguente uguaglian- 
za tra vettori: (1,0,0) = an(l,l,0) + & 2 i( 0 , 0 ,l) + 1 , 0, 1). Quindi 

(&n, &21, &3i) sono le soluzioni del seguente sistema lineare non omogeneo 
(scritto in forma matriciale): 


1 0 - 1 \ fa n \ 

10 0 a 21 

oii/ vw 



Analogamente per gli altri vettori della base canonica. Svolgendo i calcoli si 
ottiene la matrice di passaggio : 






Sostanzialmente, la difficoltà risiede tutta nel determinare le nuove coordinate 
degli elementi della vecchia base. 


Osservazione 5.5.5 Supponiamo di conoscere la matrice di passaggio dalla 
base B alla base B' . Vogliamo ricavare la matrice di passaggio dalla ba- 
se B' alla base B, cioè la matrice Mg'. Per come abbiamo definito le co- 
se, se (xi, . . . ,x n ) sono le coordinate di un vettore v rispetto alla base B e 
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(x [, . . . , x' n ) le coordinate di v rispetto alla base B r , valgono le due seguenti 
uguaglianze: 


( 


( x 1 \ 

• 

• 

• 

~ Mf, 

• 

0 

0 

vJ 


\X n ) 

( 


( ^ ^ 

• 

0 

0 

Mi' 

• 

0 

0 

\Xn) 


{<) 


In definitiva, per ogni n-upla (xi, . . . ,x n ) deve valere: 

( %i\ / ^i\ 


= Mf M| 





Ma questo implica che • M|, = I n , dove I n è la matrice identità di 
ordine n. Quindi, M| è la matrice inversa di Mg,. In particolare, le matrici 
di passaggio sono matrici invertibili. 


Esempio 5.5.6 Consideriamo lo spazio M 2 (M) delle matrici quadrate di or- 
dine 2. Questo è uno spazio vettoriale di dimensione 4 e le quattro matrici 

A\, A 2 , A$, A 4 : 




7 


costituiscono una sua base. 

In questa base, la generica matrice A = 


a 


c d 


ha coordinate (a,b,c, d). 


Consideriamo ora il sottospazio vettoriale S di M 2 (M) delle matrici sim- 
metriche, 

S = {A <E M 2 (R) | l A = A} . 


Si può verificare (esercizio) che dim S = 3 e che una sua base è data dalle 
tre matrici: 
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Completiamola ad una base di M 2 (M) aggiungendo la matrice 



e consideriamo la base ordinata Ai, E 1? E 2 , A 4 . 

Chiaramente, le nuove coordinate di Ai e di A 4 coincidono con le vecchie 
e sono, rispettivamente, (1,0, 0,0) e (0,0, 0,1). Osserviamo che: 



1 

2 


Ei H — E 2 
2 


e A 3 



Quindi le nuove coordinate di A 2 e di A 3 sono, rispettivamente : (0, |, |,0) 
e (0, 0). 

A questo punto è facile costruire la matrice di passaggio; basta scrivere 
la matrice che abbia come colonne le nuove coordinate di Ai, A 2 , A 3 , A 4 : 



fi 

0 

0 

°\ 

0 

1 

2 

1 

2 

0 

0 

1 

2 

1 

2 

0 

\0 

0 

0 

V 


In questa nuova base, la generica matrica A 




ha coordinate 


/xi\ 

( a \ 

/ a \ 

x 2 

x 3 

= M 

b 

c 


Ò+C 

2 

b—c 

2 

W 

w 

\dj 


5.6 Somma e somma diretta di sottospazi vet- 
toriali 

Proposizione 5.6.1 Siano V e W due sottospazi vettoriali di uno spazio 
vettoriale U . Allora V fi W è anch’esso un sottospazio vettoriale di U . 

Dim. Dobbiamo dimostrare che V fi W verifica le condizioni a) e b) della 
Definizione 5.2.1, sapendo che sia V che W le verificano. Siano x, y E V C\W . 
Gli elementi x + y, \x appartengono a V perché V è un sottospazio; per lo 
stesso motivo, x + y, \x G W. Allora sia x + y che \x appartengono a V fi W. 

□ 
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CAPITOLO 5. SPAZI VETTORIALI 


In particolare l’intersezione tra due sottospazi vettoriali non è mai vuota. 
Essa contiene almeno lo zero. 

Dati due sottospazi V e VE di uno stesso spazio vettoriale [/, l’unione 
V U EE non è, in generale, un sottospazio vettoriale. 

Prendiamo come esempio 


X = {(x,y) el 2 



0} e Y = {(x, y ) e M 2 


x = 0} . 


I vettori (1,0) e (0, 1) appartengono alUF ma la loro somma (1, 1) no. 
Dunque X U Y non è un sottospazio vettoriale di R 2 . 

E possibile però costruire a partire da E e VE un altro sottospazio di U 
nel modo che segue: 


Definizione 5.6.2 Si dice somma dei sottospazi V e W di U il sottospazio 


V + VE = {x eU 


x 


v + w,v £ V e w £ VE} . 


Questo è l’insieme di tutti i vettori di U che possono essere scritti come 
somma di un elemento di V più un elemento di VE. 

Possiamo dimostrare che: 

1) l’insieme V + VE è effettivamente un sottospazio vettoriale di U ; 

2) V + VE contiene come sottospazi sia V che VE ; 

3) V + VE e il più piccolo sottospazio di U che contenga sia V che VE. Ov- 
verosia, se X è un sottospazio di U tale che V C X e VE C X, allora 
necessariamente V + VE C X. 

La dimostrazione di questi fatti viene lasciata per esercizio. 

Abbiamo visto, dalla definizione, che ogni elemento di V + VE si scrive 
come somma di un elemento di V con un elemento di VE. In simboli: 


x <E V + W 3v e V e 3w e W tali che x = v + w . 

Gli elementi v e w sono univocamente determinati da x oppure possono 
esistere v' E V e w' E VE rispettivamente diversi da v e w ma tali che 
ugualmente si abbia x = v' + w'1 

Proposizione 5.6.3 Sia x un elemento di V + VE. Esiste un unico v G V 
ed un unico w G VE tali che x = v + w se e solo se V il VE = {0}. 

Dim. Supponiamo che esista z G V fi VE, z ^ 0. Allora si ha: x — v + w = 
(v + z) + (w — z) . Anche v + z^Vev + z^v. Analogamente per w e w — z. 

Viceversa sia x = v + w = v' + w' e v — v' ^ 0. Allora v + w — v' — w' = 0, 
cioè v — v' = w' — te, per cui v — v' G VE. In conclusione v — v'^VDW e 
v — v' 7^ 0. □ 


